International Journal of Electrical and Computer Engineering (IJECE) 
Vol. 12, No. 3, June 2022, pp. 2393~2404 
ISSN: 2088-8708, DOI: 10.1159 1/ijece.v1213.pp2393-2404 o 2393 


A sensorless approach for tracking control problem of 
tubular linear synchronous motor 


Nguyen Hong Quang!, Nguyen Phung Quang”, Nguyen Van Lanh? 
l Faculty of Mechanical, Electrical, Electronics Technology, Thai Nguyen University of Technology, Thai Nguyen City, Vietnam 
Institute for Control Engineering and Automation, Hanoi University of Science and Technology, Hanoi, Vietnam 
3School of Electrical and Data Engineering, Faculty of Engineering and Information Technology, University of Technology Sydney, 
Sydney, Australia 


Article Info ABSTRACT 
Article history: As well-known, linear motors are widely applied to various industrial appli- 
Received Jul 16. 2020 cations due to their abilities in providing directly straight movement without 


auxiliary mechanical transmissions. This paper addresses the sensorless con- 
trol problem of tubular linear synchronous motors, which belong to a family 
of permanent magnet linear motor. To be specific, a novel velocity observer is 
proposed to deal with an unmeasurable velocity problem, and asymptotic con- 
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1. INTRODUCTION 

During the past few years, tubular linear synchronous motors (TLSM) have been intensively applied to 
various applications including robotics, transportation, and vehicle systems because of the absence of mechan- 
ical reduction and transmission devices (gears and lead screws) permits to obtain higher precision and reduced 
dimensions in comparison with rotary motors. Recently, the merits of using TLSM have been pointed out by 
El-B], which can be listed as low cost, the durable structure, reliable operation. In addition, there has been a 
large number of researches devoted to applications of TLSM such as active vehicle suspension [4], the planar 
magnet array [5], jetting dispenser [6], and two-dimensional nanopositioning [7]. More recently, the progress 
of modern control engineering has been toward the tracking problem of TLSM, which can be noted as thrust 
optimization [B], model predictive control i), (10), and fuzzy control mi, (12). 

Tubular linear synchronous motors whose structure as shown in Figure | contains a tube (slider) with 
mounted drive magnets, and three phases winding in stator placed differently 120° of electrical angle. Without 
auxiliary reducer or transmission, TLSM is capable of operating effectively by eliminating mechanical hystere- 
sis. That however also rises the sensitiveness on the movement of slider due to frictional force, load variation, 
and non-sinusoidal flux. These unexpected forces shrink performance of the motion system both the transver- 
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sal and in the longitudinal direction. In recent years, there has been a majority of works devoted to improving 
the position tracking the performance of linear motor systems under the impact of external disturbance [13]. 
Besides, a backstepping sliding mode controller based on nonlinear disturbance observers was created by [14], 
aimed at obtaining tracking performance as well as disturbance rejection. An adaptive robust controller 
was proposed by [16], in which the dead-zone compensation technique is applied to guarantees tracking perfor- 
mance and system robust against uncertainties. Furthermore, the researchers in present an effective neural 
network learning controller for tracking the position of the linear motor. 
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Figure 1. Structure of tubular linear motor 


In recent years, sensorless control for linear permanent magnet motor have received a great deal of at- 
tention from industrial applications including: linear tubular motors [12], [18], [19], position sensorless control 
(4), (20), RI], end-effects [11], [22], [23]. Most of the researches on sensorless approach for linear perma- 
nent magnet motors has taken advantage of back electromotive (EMF) [19], [24], currents, and voltages into 
account to observe the velocity. Unfortunately, the main drawbacks of these methods are that the measure- 
ment of currents and voltage usually contain unexpected noise due to the impact of a power converter [25]. 
Moreover, at low speed or stopping operation, the value of EMF could be unreliable or vanish, which leads to 
the difficult implementation of the sensorless approach. Further, the occurrence of uncertainty in parameters 
(inductance and resistance) of the TLSM may result in noticeable estimation errors. Besides, to compensate 
impacts on control performance of the uncertainties and modeling imprecision (e.g. friction, parameter varia- 
tion, and load disturbance), the adaptive method [26], is employed to determine the required thrust force. 
The studies however have not been concerned with the sensorless control problem. It is worth noting that 
the position sensor is usually mounted on TLSM and the ordinary control scheme of TLSM is illustrated by 
Figure 2. However, obtaining velocity from the position by taking a derivative can lead to inaccurate results 
due to the noise contained in a position measurement. 

From the issues pointed out above, this paper concerns tracking control problems of TLSM subjected 
to unknown loads. By employing a novel observer-based control approach, our work aims to solve a practical 
problem that the velocity of TLSM can not be measured directly and differential calculation from the measured 
position is inaccurate. To outline, our contributions can be highlighted: 

- As mentioned above, the use of sensorless approach involves difficulties in dealing with noise measure- 
ments and electrical parametric uncertainties. To alleviate these concerns, this paper provides a novel 
observer which requires the measured position to estimate the velocity of TLSM. This observer ensures 
that the velocity error exponentially converge to zero. Moreover, the method of choosing proper param- 
eters for the observer is presented. 

- In addition, the multiple-loop control design which includes both a position-velocity controller and a 
current controller is provided to improve performance of the tracking control problem. By using the 
Lyapunov direct method, the position the position and velocity tracking errors are ensured to converge to 
small arbitrary values. 

The organization of this paper includes 5 parts as follows. The next presents a dynamic model of the 
TLSM in d—q axis and shows the main problems in sensorless control of TLSM. In section 3, a unique velocity 
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observer is introduced, and asymptotic convergence of observer errors are proofed. After that, in section 4 the 
position-velocity tracking control system for TLSM is proposed. Later, the verification of whole system is 
demonstrated by simulation results in section 5. Finally, conclusions are summed up in section 6. 
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Figure 2. Typical field-oriented control (FOC) diagram of TLSM 


2. PROBLEM STATEMENT 
As aforementioned, the TLSM has three separated wingdings a, b, c which contain three AC currents 
as: 


2 
ia) = Tsin (wt), i(t) = Tsin (wt F), it) = Tsin (we +), 


where w is electrical velocity of the motor. Notice that the electrical angular position of TLSM can be expressed 
with respect to primary position as (1) 


6(t) = “a(t), w(t) = Tufi). (1) 


By using Clarke and Park transform combine with (ip, the dynamic model of tubular linear motor is expressed 
in d — q as (2)-(5): 


Z Sof), (2) 
T = aight) — — felt) — Ill) 3) 
di Rs. 27 ; 

na = a(t) 4 ( í u(t) ig(t) + ualt), (4) 
di Rs. 27 . 2T 

at =~ (0) - (Foe) iat - (Fo) P+ Fu, © 


where iq(t), iq(t) are the stator current projected in d, q axis; the linear velocity and position of rotor is denoted 
as u(t), a(t); wp, is the flux of the permanent magnet; Rs, L stand for, respectively, stator’s resistance and 
inductance; m, T represent for the slider’s mass and pole step length. The inputs voltage of system in d — q axis 
are denoted as uq(t) and u(t). The disturbance consists two factor, the first one is the disturbance load f¢(t), 
the other one is the force generate by inductance fluctuation combines with the detent force which is 
represented as fm (t). 

To be specific, the reference position and velocity of the TLSM are denoted as x(t), v,(t). The goal 
of this research is to control both position and velocity of the TLSM, by which the actual position and velocity 
follow the desired trajectory, x,.(¢) and v,(t), with desired small errors and robust against the load variations. 
In fact, rotation motors can be easily setup velocity measurements by attaching a encoder or resolver, that of 
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linear motors however is challenging due to mount, high cost and sensitive with external factor like humidity, 
temperature, vibration. To alleviate these concerns, a novel observer is designed to estimate the velocity from 
available data of the position and current sensor. In addition to concern reluctance effects and force ripple 
reduction, ig(t) should be regulated to the reference iar (t) = 0. 


3. VELOCITY OBSERVER DESIGN 
Shortly, let us define fa(t) as the sum of disturbance i.e fa(t) = + fe(t) + + fm(t). Furthermore, the 
disturbance is assumed to be bounded such that 


| fa(t)| < Fm, | fa(t)| < dF, (6) 


where Fm and dFm are given positive constants. The following observer plays a key role in the derivation of 
our approach 


A(t) = 0(¢) + pe(x(t) — (0), r 
A 37 Wp a 
je) = ype ioe #(t)) + ysign(a(t) — #(4)), 
in which pz, py, y are real positive constants, and #(0) = 2(0). Let ay )= AUI Z(t), v(t) = v(t) — ê(t), 
then the observer errors dynamics can be obtained by the help of (2), (3) and (7) as (8) 


w(t) = —peB(t) + (t), 
v(t) = =p? (t) — ysign(Z(¢)) — falt). 
The following theorem provides a choice of observer parameters pz, pu, y by which the observer (8) are well- 


posed. 
Theorem 1: For proper positive constants a > 0, let pz, py, y satisfying 


(8) 


Pepy 0 Put 502 ipe 
> 2a oe AR ; 9 
| 0 Pr | z | 3Pr 1 a) 
T sEm 
TE — Ee — dFin > aly + Fin), (10) 


the system (8) is exponentially stable. Moreover, there exist c1, c2 > 0 such that 
|z(t)| < cae, |- PrE) +) < ee. 


Proof of Theorem 1: Denoting C1 (t) = Z(t), &2(t) = —px%(t) + (t) and C(t) = [G (t), &2(t)]T, then (8) can 
be rewritten as (11) 


Ġ& (t) = &(t), 
C2(t) = — put (t) +Ù) (11) 
= — pakolt) — Po& (t) — ysign(¢1 (¢)) — falt). 


Examine the function as (12) 


Hct) = 580 + (2+) ROHANO MOG + Fae). a 


And recalling that y > Fm, we have y|¢i(t)| > | fa(t)Gr(0)]. Then, V;(t,¢(t)) is a positive real function, 
furthermore it has Vi (t,0) = 0, Vi(t, ¢(#)) > OV¢(t) Æ 0, and Vi (t, ¢(t)) oo as ||¢(t)|| — oo. According 
to solution of (Ip, let us take the time derivative of as (13) 


HEED) = GOGO + (pe + 2) GOO HAGO) +O 
+ fa) + FEO + FaWAW 
=-F GH + faa) + FAO hE) -signa lt) fal) 


= Eee) — P(e) + a(t) (sin(a) -E +f). a 
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Using y in theorem 1, it is worth mentioning that 


We > | PE falt) + ja): (14) 
In the light of (14), results in 
z Px 2 Px Pv 2 YPx Palm 
Vit, C(t) <— Ego -Eao — (Ee Pe dn JIG). (15) 


From (10), is ensured by 
1,2 1 
wince) < 57 | PEA Be | c+ EAO 


By recalling conditions b} and {10}, it leads to Vi (t, C(t)) < —2aVj(t, C(t)). Then, using comparison lemma 
in Lemma 3.4 and the initial condition 7(0) = ¢1 (0) = 0, we have that 


Vi CE, CC) < Va(0,C(O))e-28* = S(O)?" 


Accordingly, Çı (t) and ¢2(t) exponentially converge. Then, exist c1, c2 > 0 such that ¢? (t) < ciVi(t, ¢(t)), 
G3 (t) < c3Vi(t, C(t)). Obviously, the theorem 1 is proved. 

Remark 1: With the assumption of no load disturbance (fa(t) = 0), the switching term in (7) is no longer 
needed. By excluding the switching term, the proposed observer be become the high-gain observer as in [30]. 
Hence, the observer (7) can be seen as an improvement for high-gain observer that address to handle the impact 
of disturbance. 


4. CONTROLLER DESIGN 

Using a cascade control strategy, we separate the TLSM system as presented in (2)-(5p into two sub- 
systems which are position-velocity (outer subsystem) and current subsystem (inner subsystem). It should be 
noted that the time response of inner subsystem is much faster than that of the outer subsystem. The two control 
loops are present as follows. 


4.1. Velocity-position controller 
For simplicity sake, The desired velocity and acceleration is respectively denoted as v,(t) = «,.(t), 
a,(t) = ù, (t). Further, we define these following symbols 


o= ae ex(t) = x(t) = x,(t), (16) 


ev(t) = v(t) — vr (t), Ev (t) = Ot) — v(t), 
And follow the position-velocity subsystem in 2)-@) becomes 


s(t) = ex(t), 


é,(t) = oi" — ap(t) — fall), os 


In which i; stands for the reference quadrature current which is assigned to inner control loop. Apply the 
assumption that i,(¢) simultaneously track 7}. Then, replaced i, with i} (t) in (17. The controller for outer 
loop is provided as (18) 


Theorem 2: Consider 
z% I ` 
tq z (a(t) > Reet) = Eee (18) 


If kz, ky € R44 are large enough constants, then, the outer loop is stable, and e,.(t), e,(t) converge to 
arbitrary small values. 
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Proof of Theorem 2: Using notations in (16), it lead to é,(t) = e,(t) — v(t). Then, can be rewritten as 
(19) 

s(t) = ev(t), 

v(t) = —kzex(t) — kyey(t) + kvlt) — falt). 


On studying the control performance and stability of the closed loop system, a Lyapunov candidate function is 
applied as 


(19) 


V2(t) = 
By using (19}, it establishes that 
Va(t) =(ke + kv )€x(t)ev(t) + ev (t)éu(t) + eb (t) + ex (t)év(t) + Vi(t) 
= — kee2(t) — (ky — 1)e3(t) + koB(t)(ea(t) + ev(t)) — fault) (colt) teo(t)) V4). 20) 
From (15), it is clear that V; (t) < — 22 (Ps1 (t) + Ga(t))? = %0 (t). By applying the inequalities as 


IT) (es (t) + ex (t))| +ee(t)+ (E = )*0, 


(ke + kode) + 5e2(t) + ealt) + Vall). 


cone 


| fa(t) (ex(t) + ev(#))| < iele 2(t) + e?(t)), 
where €x, €v, €f > 0, from we obtain that 


Va(t) < — (ka — €x — ef) (t) — (ky — &» — € — Lee (2) 


7 1 TV F? (21) 
P ve (t) it m . 
8 de, de, Qe Fr 


Accordingly, by choosing 
ky — €, — ef = 1, 


ky — ev —ef —1=1, (22) 


1 _ 


2px — Ez TS 0. 


Then V2(t) < 0 for all (ex(t), ex (t)) ¢ E where 


F2 

E£ flene) ER: ete <— ar (23) 
2ef 

It implies that the tracking errors (e(t), ez(t)) enter £ in finite time due to V2(t) < 0. By choosing ep large 

enough, the tracking errors can converges to arbitrary small values. Intuitively, the theorem 2 is proved. 


4.2. Current controller 

From the fact that dynamics of the current loop is always much faster than that of outer loop, the 
reference i; can be assumed to be unvarying in inner-loop control process. Additionally, the inconstancy 
in inductance cause by end-effect phenomenon can be neglected. Continuously, the following notations are 


provided 
Cig(t) = tq(t) — 9, eia(t) = ta(t) — iå 
In what follows, let us establish a modified PI controller which cooperates with the velocity observer in 


section 3. The current loop can be consider as two parallel current system and therefore can be controlled 
by two PI-like controller as (24), (25) 


ualt) = Rai — hacialt) = kia f eaat - ig OOO, (24) 
g(t) = Rait — kgeig(t) — kiq | eig(E)dé + Gaz + He | a(t). (25) 
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In which ka, Kp, kia, kiq are positive constants. From that, the inner closed-loop is derived 


ċia(t) = at Ae) i f eial€)dé + TOT), (26) 
eit) =- eult -e f eua- (Zi + Tee). 27 


Using the same method which presented in subsection 4.1, by applying the Lyapunov candidate function as in 
(28) 


t 2 l t 2 
Vile) = Vil AO O+ E f eaa) ieh f eoa). es) 


It is clear that the tracking errors of both quadrature and direct current converge to zero. In concluded, the 
control scheme whole system is describe in Figure 3. 


Current 
Controller 


Figure 3. The proposed control diagram for the TLSM 


5. SIMULATION RESULT AND ANALYSIS 

TLSM is used in this simulation has the parameters as follow: Rs = 10.3 (Ohm); m = 0.171 (kg); 
L = 1.4 (mH); np = 1;7 = 10 (mm); Yp = 0.035 (Wb). This motor is mounted a position sensor in which 
output data are sin(@) and cos(@). The control performances of the sensorless approach are verified by two 
difference simulation scenarios. On one hand, the first scenario verifies the effectiveness of the control system 
in the case of no measurement noise affect the TLSM. In the other hand, the other test is done with the impact 
of measurement noise to the position feedback signal. In both scenarios, the using controller, observer and 
disturbance load are the same. 

The simulation of the observer based control for TLSM is ranon MATLAB/Simulink with the chosen 
sampling time is 10~°(s). The load’s disturbance applied to TLSM given by: fe(t) = 3 + 48 sin(20¢) + 
32 sin(60t) + 32 sin(100¢). It is worth mentioning that, in practical case, fm (t) is very small in compare with 
fe(t). Therefore, we take Fm = 60,dF,, = 2000. From (9). {10}, observer’s parameters are given by 
Px = 10°, p, = 2.104, y = 100, a = 30. Following (22}, the parameters of controller block are chosen as 
Kp = 105, k, = 2.108, ka = Kp = 10, kia = kiq = 104. 


5.1. None measurement noise case 
As mention above, the feedback position from sensor is assumed perfectly accurate, the initial errors 
of observed position and velocity are chosen as %(0) = 0, 0(0) = 0.1. As illustrated in Figures 4(a) and 4(b), 
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the actual position and actual velocity of TLSM follow the desired trajectories. Also, there are variations in 
tracking errors in Figures 4(c) and 4(d) during the response period of the observer. Values of the current and 
voltage of the q-axis are shown in Figures 5(a) and 5(b). Although the error still converges to zero in less 
than 0.1 s Figure 5(c), the TLSM is affected by the disturbance load as in Figure 5(d). During this interval, 
the motor stay still, which is the advantage of propose velocity observer compare with the other EMF-based 
techniques. Accordingly, theorem | is verified. Furthermore, under the disturbance load, the position tracking 
is still maintained. Figure 5(a) depicts that the actual quadrature current i,(t) follows the desired signal i 
With the high precision tracking and quick response time, the simulated results confirm the preformance of the 
propose control system for TLSM. 
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Figure 4. Position and velocity behavior of the TLSM, (a) position, (b) velocity (c) position error, 
and (d) velocity error 


5.2. Presence of measurement noise case 
In this case, opposition from the first case, The data from position sensor is assumed to contain mea- 
surement noise as (29) 


Lmeas(t) = x(t) + n(t), (29) 


in which £meas(t) denotes the measured signal from the position sensor, and n(t) is a noise measurement 
which is normally modelled by a white process. Difference from the first case, the initial error of observed 
position is not equal to zero due to the measurement noise. Therefore, the initial errors of observed posi- 
tion and velocity are chosen as Z(0) = 5.1073, 0(0) = 0.1. With the measurement signal as depicted in 
Figure 6(a), the widely used method which consist low-pass filter combine with derivatives can not obtain the 
accurate velocity. The voltage value is illustrated in Figure 6(b). Overall, the proposed controller still outweigh 
in position tracking control, as depicted in Figures 6(c) and 6(d). Current and velocity values are shown in 
Figures 7(a) and 7(b). As illustrated in Figures 7(c) and 7(d), the errors between actual and observed values 
converge in approximately 0.1 s and no greater than 0.05 (m/s) with velocity and 0.002 (m) with position 


errors. 
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Figure 5. Time response of the TLSM without measurement noise, (a) current i4, (b) voltage uq, (c) velocity 


observer error, and (d) disturbance load 


0.05 ; f 4) i 
x(t) v(t) 
SRR Xr(t) ~ 
-~ ” 
E E 2| 
5 z 
2 8 
a 2° 
-0.05 2 
02 04 06 08 1 0 02 04 06 08 1 
Time (s) Time (s) 
a, (b) 
10 %10" . i - 
ex(t) T 0.05} evtt) | 
E E 
5 5 5 o 
= - 
v o 
S 2 0.05 
= Li 
3° $ 
> 
E -0.1 
-5 L — a a = 
0 02 04 06 08 1 o 02 04 06 08 1 
Time (s) Time (s) 
(c) (d) 


Figure 6. Position and velocity behavior of the TLSM under noise measurement, (a) position, (b) velocity, 


(c) position error, and (d) velocity error 
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Figure 7. Time response of the TLSM under measurement noise, (a) current 7,, (b) voltage uq, (c) velocity 
observer error, and (d) position observer error 


Remark 2: Under the impact of measurement noise n(t), the traditional method which obtain the velocity 
via a differential calculation of measured position can not estimate the actual velocity. It should be noted 
that employing low-pass filters on measured signal seems to be ineffective because the measurement noise 
often is not considered a deterministic signal. Additionally, the use of Kalman filter [32], in this case 
possibly results in the large velocity estimation error due to the effects of unknown disturbance load in velocity 
dynamics. Further, at low speed or stopping operation, the value of EMF could be unreliable or vanish, the 
EMF approach [19], leads to the difficult and inaccurate velocity estimation. From the above analyses of 
previous approaches, the proposed method shows advantages in the velocity estimation in the occurrence of 
measurement noise n(t). 


6. CONCLUSION 

This note has provided a novel technical solution for the sensorless tracking control problem of TLSM 
under the lack of velocity sensors and unknown disturbance loads. The main contributions of our method 
have based on the proposed velocity-observer, which ensures asymptotic the convergence of observer errors. 
By cooperating with the observer, the position-velocity tracking controller and current controllers have been 
constructed by using Lyapunov direct method. These controllers have ensured that the position and velocity 
error converges to arbitrarily small values by choosing properly control parameters. In later work, the current 
sensorless control will be taken into account with no further sensor requirement in control TLSM. 
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